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Abstract 

. We derive determinant representations and nonlinear differential equations for 

I the scaled 2-point functions of the 2D Ising model on the cylinder. These equations 

CS| ' generalize well-known results on the infinite lattice (Painleve III equation and the 

^ ■ equation for the r- function of Painleve V) . 

> ■ 

m ■ 

"-^ '. 1 Introduction 
o . 

oo ■ 

O ' One of the most beautiful results of the modern mathematical physics is the connection 
^ ■ between correlations in the 2D Ising model and Painleve functions. Quantum field theory 
, and condensed matter physics usually deal with Dyson equations for the Green functions. 
^— > ! These equations, however, can not be reduced to closed form, i. e. to find a propagator 
Q^l from Dyson equations one must first perturbatively calculate the vertex. The Ising model 
provides a remarkable example of closed differential equations for the correlation functions 
and, moreover, these equations prove to be integrable. 
•rH , Scaled 2-point function of the 2D Ising model on the infinite lattice 

^' 

(a(0,0)a(r,,r,))(-)=er_(t), s>l (1.1) 

(a(0,0)a(r,,r,))W=er+(t), s<l (1.2) 
can be expressed through the functions Tj-{t) of the single parameter 



Is- 11 



t=^^x/rl + rl (1.3) 



where 

s = sinh2/C, ^=|l-s-^|^ (1.4) 

and /C is the Ising coupling constant. 
As it was shown in [H [2] , the function 



1 



satisfies the equation 
and the asymptotic boundary conditions 



v" = -ivf-lv' + v'-V-\ (1.6) 

7] U 



V{0) ~ 1 - -Ko{2e) (1.7) 
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as 6 ^ oo. The equation (1.6) is a particular case of Painleve III equation. Furthermore, 
if we define 

then in both cases [3] 

{tc'f = 4 {tC - 0' - 4 icf {tC - C) + (CO' ■ (1-9) 

This is the equation for the r-function of Painleve V. The most interesting feature of (1.9) 
is that it contains correlation functions t± themselves instead of their combination (1.5). 

There are several similar results for the correlations in one- dimensional quantum mod- 
els: impenetrable Bose gas [1], XY and XXZ spin chains O [6], sine-Gordon field model 
[7j. In two dimensions, there were no examples of this kind other than spin-spin corre- 
lation functions of the Ising model on the infinite plane. Quite recently [U [9] the form 
factor representation for the correlation function was calculated for the Ising model on 
a cylinder. This gave an opportunity to compute the susceptibility and to investigate 
its singularity structure in the complex temperature plane in the case of a finite-size lat- 
tice. In the present paper, we continue the study of correlations in the Ising model on a 
cylinder. 

In this case, the rotational invariance is broken even in the scaling limit and one has 
three scaling parameters: 

where N is the number of sites on the base of the cylinder. Correlation function in the 
scaling limit is given by 

(a(0,0)a(r.,r,))(-) =eeT(/3)e-l^l/^(^)r-(^,l/,/3), ^ > 1, 

(a(0,0)a(r.,r,))W =eeT(/3)e-l^l/^(^)r+(x,i/,/3), s < 1. 

The notation will be explained in the next sections. Before summarize the main results 
of the paper, we define complex variables 

\x\ -iy _ \x\ +iy 
z = , z = (l-lO) 



and the function 

= In ^~ ^+ . (1.11) 

r_ - r I 
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It will be shown that Lp satisfies the sinh-Gordon equation 

V9^2 = i sinh 2v9. (1.12) 

Besides that, we will derive the equations for t± separately. If we define 

M = In r±, (1-13) 

then u satisfies the equation 

in both cases. This is quite a fascinating result. In fact, we will find a general A^-solitonic 
solution of (1.14) and the Ising scaled correlation function is an infinite-solitonic one. 
As far as the author is aware, the equation (1.14) is not one of common knowledge. 
Nevertheless, it seems to be integrable and complete investigation of it (construction of 
the Lax representation, hamiltonian structure, etc.) is a challenging task. 

This paper is organized as follows. The scaling limit of the formulas of [H [9j and 
determinant representations for the correlation functions on a cylinder are obtained in 
Section 2. On the infinite lattice correlation functions can be expressed through Fredholm 
determinants of integral operators and on the cylinder — through the determinants of 
some infinite dimensional matrices. The derivation of the differential equations (1.12), 
(1.14) and a brief rederivation of Painleve III (1.6) and Painleve V (1.9) equations are 
given in Section 3. 



2 Scaling limit and determinant representations 
for the correlation functions 

The most convenient representations for the investigation of correlation functions on the 
infinite lattice are the so-called form factor expansions pJi \TT\ : 

00 

(a(0,0)or(r,,r,))(-) =e^(?2n, s > 1 (2.1) 

n=0 

00 

(a(0,0)a(r,,r,))W=eX;^?2n+i, s<l (2.2) 

n=0 

Iff Q-\rx\'yilj)+iryqj 

9- = ^^J ^^^■■■^^-^^^m sinh7(g,) ' = ^'-'^ 

— TT —TV 

l<i<j<n 2 
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where the function 7(g) is defined by the equation 

cosh 7(g) = s + s^^ — cosg 

and condition 7(g) > 0. 

It proved [8l |9] that the correlation function on the cyhnder with sites on the base 
may be written in a similar form, namely 

[iV/2] 

(a(0,0)a(r.,r,))(-) = ■ ^t) e'l^^l/^ Yl 92n, s > 1 (2.5) 

n=0 
[{JV-l)/2] 

(a(0,OMr,,r,))W = (^CT)e-l^'-l/^ J2 92n+i, s>l (2.6) 

n=0 

p-'^/A___(M " p-krcl7(gj)+*»-y'?j-'^(gj) 

[g] j = l '^^^^ 

The main differences between (2.1)-(2.3) and (2.5)-(2.7) are the substitution of the 
integration by the summation over discrete bosonic spectrum of quasimomentum {qj = 
k, k = 0, 1, . . . , — 1) and the appearance of cylindrical parameters ^t, A and z^(g) 



2n 
N 



f dp dQl' {p)l' (q) 



^^^^ 27r2 J sinh(A^7(p))sinh(A^7(g)) 



sin ((p + g)/2) 



(2.8) 



sin((p-g)/2) 


TV 

A-^ = ^ j dp In coth (A^7(p)/2) , (2.9) 


TT 

1 f dx) (cost) — G^'^^^'^^ 

u(q) = - ^ \ ^lncoth(A^7(p)/2). (2.10) 

TT J cosh 7(g) — cosp 



The finite number of terms in the sums (2.5) and (2.6) provides an opportunity of inde- 
pendent check of the formulas [H [9] for small N. 

The scaling limit of the Ising model is of interest mainly due to its applications in 
quantum field theory. This limit implies that 

7(0) — > 0, — > 00, Vy ^ 00, N 00, 

7(0)ra; — > X = const, 'y{0)ry y = const, 7(0)A^ [3 = const. 
In this case |12j the formulas (2.5)-(2.7) reduce to 

00 

(a(0,0)a(r.,r,))(-) = eeT(/?)e-l^'l/^(^) Y^^^n. s > 1, (2.11) 

n=0 
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(a(0,0y(r,,r,))(+)=eeT(/3)e- 



\x\/m 



^5'2n+l, S < 1, 
n=0 



n 



n -|x'U/l+(^)'+iy^-i^(«,-,/3) 



1 + 



2-kL 



13 /3 



l<i<j<n /| _|_ j 27r/j 



/3 ; 



The summation over [/] 



oo oo 



E = E ••• E 

[i] = — oo il= — oo in = — oo 



(2.12) 



-, ^0 = 1, (2.13) 



(2.14) 



is taken over integer values of Ij, and the quantities ^t(/3), A(/3), I'll, (3) are defined by 
the formulas 



27r2 




p dp qdq In 



/ / Vp' + 1 V?' + 1 sinh + l) sinh (^p^q^ + l) 



TT 



A-i(/3) = - / dp Incoth 



TT 



1 + 



2itI 



In coth 



Let us denote 



27r/i 



Mi = 



27rZi 



+ 1 1 + 



27r;, 



It is easily verified that 



n 

l<i<j<n 



— 

Ui + Uj 



The right hand side of the last formula can be expressed through the Vandermond deter- 
minant of a matrix with dimension n x n: 



detr 



Ui + Uj 



fii^) n 

i=l \ l<i<j<n 



Ui — Uj 
Ui + Uj 
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If we make use of this fact and denote 



Ei = Ei{(3) 



1 + 



( 2-Kli \ ' 

\ ) 



(2.15) 



then for Qn we obtain 



— —TTT det„ 



EiEj exp 



\x\—iy Uj+Uj \x\+iy Uj ^ +Uj ^ 

' 2 2 2 2 



(2.16) 



Let us consider an arbitrary square matrix K with the elements K^n and construct 
the determinant 

K 



1 + 



P 



1 + - TiK + At 
p 2/32 



+ ... 



K K 

^ ^ mm »^ mn 



+ 



as the expansion over ^. After comparing these expressions with (2.16), it is easy to see 
that 



n=Q 

00 



n=0 



1--A- 



where the elements of the infinite-dimensional matrix K are obtained from 

E1E2 exp 



2 2 2 2 



Ui + U2 



(2.17) 



and li, I2 — —00, . . . , 0, . . . , 00. 

Let us summarize: correlation functions of the Ising model on the cylinder in the 
scaling limit 



(<7(0, 0)<7(r,, r,))H = eeT(/3)e-l^l/^(^) r_(x, y, /3), T < T,, 

(a(0,0)a(r.,r,))(+) =aT(/5)e-l^l/^(^)r+(x,y,/5), T > T„ 
can be expressed through the determinants of infinite-dimensional matrices: 



1 + ii^ 



(2.18) 

(2.19) 

(2.20) 
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r_ - r+ 



(2.2i; 



where K is defined by (2.17). 

On the infinite lattice K reduces to integral operator with the kernel 



K{u, v) 



\x\—iy u + v \x\+iy u +v ^ 
1 e 2 2 2 2 



TT 



U + V 



(2.22) 



acting on L^(0, cxd), and matrix determinants in the right hand side of (2.20) and (2.21) 
become Fredholm determinants. 



3 Nonlinear differential equations 

In recent years numerous connections between Fredholm determinants of certain integral 
operators and differential equations become apparent [13l [HI [15] . We will formulate here 
only necessary results. 

Let us consider integral operator K with the kernel 

E(x)E(y) , , , 

K(x,y) = ^ ' E(x) = e(x)e'^+'—. 
X + y 

It was shown in ^3] that the functions 



ip = In 
ip = In 

satisfy the differential equations 



1 + XK 


-In 


1 - \K 


1 + XK 


+ In 


1 - XK 



dzz^ = - sinh 2(p, 



dzz^ 



1 — cosh 2ip 



To derive them, one has to define 



1-^2 

K 



where Ei{x) = x'^E{x), and to prove the following recursion relations 



(3.1) 
(3.2) 

(3.3) 
(3.4) 

(3.5) 
(3.6) 



(3.7) 
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and differentiation formulas 



dzUp^q = \ + Vp^_iUq^_i) + \ + Mp,g-l) , 

dzVp^q = I (Mp,0«g,0 + ^p,0%o) + | (^p+l.g + Vp^q+l) , 
dzVp,q = \ {Up^-lUq^-l + fp,_lt;g,_l) + ^ {Vp-l^q + Vp^q^i) , 



(3.9) 



(3.10) 



dz'f = Mo,o, ^zV^ = u-i-i, (3.11) 

c}^^ = -vo,o, dzi^ = -1. (3.12) 

Two additional equations 

dzzii = - {dzvf 5 (3-13) 

d-zi; = -{dz^)\ (3.14) 



were not observed in [T^. To derive, for example, first of them, one has to differentiate 
(3.12) with respect to z and to use the recursion relation (3.8). 

It can be shown that the recursion relations (3.7), (3.8), differentiation formulas (3.9)- 
(3.12) and differential equations (3.3), (3.4), (3.13), (3.14) remain valid even when is a 
matrix (finite- or infinite-dimensional) with the elements 



CmCn e 2 



mn 



The matrix (2.17) belongs to this class, so after denoting 

¥^ = ln^— ^ = In (r!-r2), (3.15) 



r_ — r 4 



2 ' 2 
we can write down a system of equations for the correlation functions 



r_ = e 2 cosh — , (3-17) 

r+ = e2sinh^, (3.18) 

a^^v^ = ^ sinh 2v?, (3.19) 
„ , 1 — cosh 2(/3 

dz-zil^ = ^ (3.20) 

dzzil^ = - {dzif f , (3.21) 

d^z'ip = - {d^^f . (3.22) 
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On the infinite lattice in the scahng hmit correlation function and, therefore, the functions 
(/? and become invariant under rotations in the (x, y) plane. In this case, one can define 

9 = 2 (zz)^ , rj(9) = e~'^(^^'^^) and then sinh-Gordon equation (3.19) reduces to Painleve III 
equation (1.6) for the function r){9). 

Before proceed with the differential equations, let us discuss the boundary conditions 
for (3.19). It should be noted that the asymptotic conditions at \x\ — > oo in the case of the 
cylinder do not determine the solution uniquely. Indeed, the asymptotics of the functions 
T±{x, y) is determined by the first terms of the corresponding form factor expansions, and 
moreover, in the expression for gi (2.13) we should take I = (this is not the case for the 
infinite lattice, where the sum over I transforms into an integral and gives the modified 
Bessel function). Therefore, we obtain 

r+^ie-N-^(°'^), r_^l, 

as \x\ —>■ oo. The corresponding asymptotics of ip 

does not depend on y-coordinate. If we assume that the last condition defines the solution 
uniquely, it would not depend on y for all values of x and obviously, this is not true. 

If it were possible to obtain an additional differential equation, containing the deriva- 
tives with respect to ^ (this quantity corresponds to QFT temperature), one would be 
able to use the solution of Painleve III equation (1.6) with the boundary conditions (1.7) 
as the initial condition at ^ = 0. However, the derivation of such equation remains an 
open problem. This problem is rather difficult to solve, mainly due to the complexity of 
the function u{l,P). 

Therefore, we have no way other than the classical one — to determine, at least 
numerically, the function (fi{x, y) and its derivative on the whole line {y = const) or 
on the circle {x = const). Such boundary conditions, combined with the periodicity in 
y-coordinate, will completely determine the solution of sinh-Gordon equation (3.19). 

Now we can turn to the derivation of the equations for the functions r±. If one denotes 

= In r_ , 

then by means of usual differentiation and using the relations (3.17)-(3.22) it is straight- 
forward to obtain the following identities: 

cosh^-1 f dz(p \ , . 

cosh 1 f _djip\ 



cosh (fi — 1 
u^-^ = 



-1 + 



sinh J \ sinh (p 



(3.25) 
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(dz^ \ ( \ 

■ , ] +Uzz\ . , , (3.26) 
smn (f J \ smn (p j 

(dz^ \ ( dz'p \ 

■ , \+Uz-z\ . , , (3.27) 
smn (p J \sinn (p J 

( \ ^ f , 2[uzzUzz-ul^] 
Uzzzz = Uzzz ■ , + ■ , - 2ii^^ H — . (3.28) 

\?>mhif J \smhLp J cosn(y9 — 1 

The last four relations can be treated as a system of algebraic equations for three variables 
cosh^y-i ^ ^_d^_^^ (^ sinhy ) " equatiou for Correlation fuuctiou 



It, 



(it^,^ + 2ul-, - Uzz) = Uzz-zU.-^z - Uz.U^z (3.29) 



appears as the consistency condition for this system. In the case m = lnr+ all the calcu- 
lations can be performed in analogous fashion, and, though the equations (3.23)-(3.28) 
change their form, the equation (3.29) for the correlation function is the same. Since 
(3.29) is obtained from (integrable) sinh-Gordon equation (3.19) and three linear equa- 
tions (3.20)-(3.22), it is very hkely integrable. More detailed investigation of (3.29) is 
beyond the scope of this paper. 

On the infinite lattice the problem of derivation of differential equations for the cor- 
relation functions is much simpler due to rotational invariance. The functions t±, ip, ip, 
u depend on the single variable t — 2 (zz)^. Instead of (3.19) one has 

(fi" + -(fi' — sinh (f cosh (p, 

t 

and instead of (3.23), (3.25), (3.26) 

,, 1 , ((p" - 7 ip') sinh (p - ((/?') ^ cosh (p 
u u — 



t 2(cosh ip + l) 

-1 + 



„ 1 , cosh (p — 1 
u" + -u' ^ ^ 



sinh ip 



u"+'-u') =2u"(^ 
t J \smh(/7^ 

correspondingly. When we again regard the last four relations as a system for three 
variables (p, ip', ip" and if we denote 

C = tu\ 

we will obtain well-known result 

{tCf = 4 {tC - 0' - 4 (Cf {tC - C) + (CO' 

— the equation for the r-function of the Painleve V equation, — as the consistency 
condition for this system. Analogous calculations lead to the same equation for r+. 
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